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The existing 210Ph , mass balance models for the assessment of cultivated soil erosion are based on an as- 
sumption that 7!°Pb,, is quite evenly mixed within the plough layer. However, the amount of 7!°Pb,, distributed 
in the soils below the plough depth, like a downward tail in the lower part of the 7!°Pb,, profile, has been largely 
ignored. In fact, after the initial cultivation of undisturbed soils, 7!°Pb,, will diffuse downward from plough 
layer to the plough pan layer due to the concentration gradient. Assuming 7!°Pb,, inventory is constant, the 
depth distribution in the two layers of the cultivated soils will achieve a steady state after continuous cultivation 
for 10.37 years, when *!°Pb,, is evenly distributed in the soils of the plough layer with an exponential concen- 
tration decline with depth in the soils of the plough pan layer, and the 7!°Pb,, concentration at any depth will be 
invariable with time. The work reported in this paper attempts to explain the formation of the 7!°Pb,, tail in the 
soil profile below the plough depth by theoretical derivation of the 7!°Pb,, depth distribution process in the two 
layers of the cultivated soils, propose a 7!°Pb,, mass balance model considering 7!°Pb,, diffusion based on the 
existing model, and discuss the influence of the 7!°Pb,, tail to the existing model. 
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I. INTRODUCTION 


Because of temporal and spatial limitations associated with 
the application of '°’Cs technique to estimating rates of soil 
erosion and sediment redistribution [1-4], the use of 7!°Pb,, 
measurement has attracted increasing attention, as an alter- 
native approach [5-9]. Unlike '°’Cs, which is an artificial 
radionuclide that is released into the environment as a result 
of atmospheric test of nuclear weapons, *!°Pb is a natural 
product of the 7°8U decay series, derived from the decay of 
gaseous ?”’Rn, the daughter of 77°Ra [10, 11]. ??°Ra is found 
naturally in most soils and rocks and the generated 7!°Pb is in 
equilibrium with its parent. A small quantity of ??Rn diffuses 
upward from the soil, introduces 7!°Pb into the atmosphere 
and provides an input of this radionuclide to surface soils and 
sediments, which is not in equilibrium with its parent ?7°Ra. 
This fallout component is termed unsupported or excess *!°Pb 
(7!°Pb,,.) [12, 13]. Like '°7Cs, 7!°Pb,, is shown to have a 
strong affinity for soil and sediment particles. Upon reach- 
ing the soil surface, 7!°Pb,, is quickly and strongly adsorbed 
by the surface soil, and the subsequent redistribution within 
the landscape will reflect the movement of soil and sediment 
particles associated with soil erosion and sediment transport 
processes [14]. Unlike the time-dependent fallout of !37Cs, 
the atmospheric fallout of 7!°Pb,, is essentially constant for 
its natural origin [15], and the 7!°Pb,, inventory in soils can 
be invariable because the subsequent fallout provides an in- 
put to surface soils that is in equilibrium with the radioactive 
decay of the existing 7!°Pb,, in the soils. 
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In the reported works, the profiles of soils from undis- 
turbed land show that the maximum concentrations of 7!°Pb,, 
occur at the surface horizon and decline exponentially with 
depth (Fig. l(a), 1(c)), while in the profiles of cultivated 
soils the 7!°Pb,, concentration is almost uniform through- 
out the plough layer as a result of mixing caused by tillage 
(Fig. 1(b), 1(d)). For the cultivated soils, it is also obvious 
that an amount of 7!°Pb,, distributed in the soils below the 
plough depth, like a downward tail in the lower part of the 
profile (Fig. 1(b), 1(d)). The 210P tail in the profile demon- 
strates that the maximum 7!°Pb,, concentrations occur at the 
bottom of plough layer and decline exponentially with depth. 
However, the amount of 7!°Pb,.. distributed below the plough 
depth has been largely ignored to date, suggesting that 7!°Pb,, 
is quite evenly distributed and restricted to the plough layer. 
The existing *!°Pb,, mass balance models for the assessment 
of cultivated soil losses have similarly not taken into account 
the amount of 7!°Pb,, distributed below the plough depth. 

Walling and He [16] developed a mass balance model for 
estimating erosion rates on cultivated soils from 7!°Pb,,. mea- 
surement. Based on this mass balance model, and ignor- 
ing the freshly deposited *!°Pb,, removed by erosion and 
influence of the particle size correction factor, Zhang et 
al. proposed a simplified 7!°Pb,, mass balance model in 
the steady state for estimating erosion rates on cultivated 
soils [18], which can be expressed as: 


AH 


H )A=0, (1) 


I— (A+ 
where A is the 7!°Pb,, inventory (mBq/cm7), A is the decay 
constant of 210Pp (/yr), AH is the erosion rate (cm/yr) and 
H is the plough depth (cm). 

The 7!°Pb,.. loss proportion (J) relative to the local refer- 
ence inventory (Aref, in mBq/cm”) for an eroding site can be 
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Fig. 1. Typical 7!°Pb,, profiles in undisturbed and cultivated soils. 
(a), (b), Moorlake catchment, UK [16]; (c), (d), Sichuan Hilly Basin, 
China [17]. 


expressed as: 
l= [(Avep — A)/Arer] x 100%. (2) 


Assuming the *!°Pb,, deposition flux is constant, the ero- 
sion rate can be calculated as Eq. (3) according to the relation- 
ship between the deposition flux and the reference inventory 
d= AAref)- 


APs =, (3) 


So, the existing mass balance models (existing model for 
short) are based on the assumption that 7!°Pb,, is quite evenly 
distributed within the plough depth. Nevertheless, neglect- 
ing 7!°Pb,,. distributed below the plough depth, more or less, 
would consequentially influence the application of existing 
model to soil erosion estimation. So, it is needed to know the 
influence level, and whether or not the existing model is ap- 
plicable to soil erosion estimation as usual. In this paper, we 
attempt to explain the formation of 7!°Pb,, tail in the lower 
part of profile by theoretical derivation of the 7!°Pb,, depth 
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distribution process in the plough and plough pan layers of 
the cultivated soils, propose a 7!°Pb,,. mass balance model in 
consideration of the radionuclide diffusion (revision model 
for short), and discuss the influence of the *!°Pb,,, distributed 
below the plough depth on the application of existing model 
to erosion estimation. 


I. 7°PB,x DISTRIBUTION PROCESSES IN CULTIVATED 
SOILS 


A. Mechanism and process analysis 


In fact, *!°Pb,, is quite evenly mixed by tillage within 
the plough layer, and there is little *!°Pb,, distributed below 
the plough depth at the beginning of cultivation, because the 
plough depth is usually no less than the measured depth range 
of 7!°Pb,.. distribution in the undisturbed soils before tillage. 
However, this initial depth distribution state is unsteady, and 
the radionuclide moves downward from the plough layer to 
the plough pan layer. Radionuclide downward movements in 
soil include diffusion and migration processes [19]. 7!°Pb,, 
diffusion in soil is the process of radionuclide movement 
that is caused by the concentration gradient in the ionic or 
molecular form, which means that the radionuclide is only 
able to diffuse from a high concentration to a low concen- 
tration. In 7!°Pb,, migration in soil, its downward move- 
ment is caused by transportation via media, such as adsorbed 
solid particles and dissolved fluids, which means that the mi- 
gration is independent of the concentration gradient, and the 
radionuclide is able to move either from a high concentra- 
tion to a low concentration or vice versa. For a soil profile 
in which the maximum radionuclide occurs at the top hori- 
zon and declines monotonically with depth, diffusion can be 
considered the predominant radionuclide movement process 
over the profile, while the radionuclide migration processes 
can be neglected [20-24]. Thus, it can be assumed that the 
diffusion process dominates the formation of the *!°Pb,, tail 
in the profile of the plough pan layer. 

At the beginning of cultivation, 7!°Pb., is quite evenly dis- 
tributed within the plough layer; after the initial cultivation, 
the *!°Pb,, diffuses from the plough layer to the plough pan 
layer due to the concentration gradient; the 7!°Pb,,. inventory 
of the plough layer decreases while that of plough pan layer 
increases; and the 7!°Pb,, tail in the profile of the cultivated 
soils progressively comes into being closely below the inter- 
face between the two layers (Fig. 1(b), 1(d)). Assuming the 
710Pp,., inventory is constant, the *!°Pb,, concentration at any 
soil depth will be invariable after a period of diffusion, and 
the 7!°Pb.,. depth distribution in the two layers of cultivated 
soils will achieve a steady state. 


B. The depth distribution in the unsteady state 


Assuming no soil erosion or deposition occurs, when 
210Ph is evenly mixed within the plough layer by the ini- 
tial cultivation of undisturbed soils. Subsequently, the 7!°Pb,, 
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starts moving from the plough layer to the plough pan layer. 
Taking the diffusion process into account, the 7!°Pb,, move- 
ment in the soils of plough pan layer follows Fick’s second 
law of diffusion [25], which is: 


D ð (OC p(z, t) OCpp(Z, t) 
Oz Oz Ot 


ACpp(z,t) =0, (4) 


where D is the effective diffusion coefficient (cm?/yr), Cpp is 
the 7!°Pb,, concentration in the soils of the plough pan layer 
(mBq/cm*), z is the depth downward from the interface be- 
tween the plough and plough pan layers (cm), and ¢ is the time 
elapsed since 7!°Pb.. was well mixed within the plough layer 
(yr). When z — 0, the 7!°Pb,, concentration of the interface 
soils is equal to the value of the 7!°Pb,, concentration at any 
depth of the plough layer. 

By diffusion, the total *!°Pb,, (A, in mBq/cm?) is redis- 
tributed between the plough and plough pan layers, which 
means: 


A = Ap + App, (5) 


where A, and App are the 7!°Pb,, inventory of the plough 
layer and plough pan layer, respectively. 

The 7!°Pb,, inventory in the soils of the plough layer can 
be calculated as: 


Ap = HCy. (6) 


Assuming the plough pan layer can be treated as a semi- 
infinite homogeneous porous medium, the process of 7!°Pb,, 
diffusion from the plough layer to the plough pan layer can 
be considered as the infinite thin layer radionuclide resource 
diffusing to a one-dimensional semi-infinite medium. Thus, 
according to the Gauss solution of the Fick’s second law [26], 
by combining Eqs. (4), (5) and (6), the sea 3s concentration 
variation with time in the soils of the plough layer can be 
expressed as: 


Cp(t) = A/[H + (nDt)'/?]. (7) 


And the variation of the 7!°Pb,, concentration with time 
and depth in the soils of the plough pan layer can be expressed 
as: 


A 22 
— g iD, (8) 
H + vnDt 
By integrating the two sides of Eq. (8) over depth z, the 
proportion of the 7!°Pb,, inventory of the plough pan layer to 


the total 7!°Pb,, inventory in the cultivated soil, namely the 
210Pb... diffusion proportion, can be expressed as: 


App(t) vaDt 
A H + vnDt 
The above theoretical derivation, namely Eqs. (4)-(9), can 


also apply to the '°7Cs depth distribution process in cultivated 
soils. However, unlike '3’Cs, the 7!°Pb.,. depth distribution in 


Cw (z, t)= 


x 100%. (9) 
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the two layers of the cultivated soils finally achieves a steady 
state after a period of diffusion, assuming the *!°Pb,, deposi- 
tion is constant over time. Before the steady state, the 7!°Pb,, 
concentration at any depth of the two layers varies with time. 
If the effective diffusion coefficient D is known (which can 
be obtained by experiments), Eq. (7) can be used to simulate 
the 7!Pb., concentration variation with time in the soils of the 
plough layer, and Eq. (8) can be used to simulate the 7!°Pb.., 
concentration variation with time and depth in the plough pan 
layer. 


C. The depth distribution in the steady state 


After achieving a steady state, the !°Pb,, concentration at 
any depth of the cultivated soils is invariable with time, viz. 
OCpp(z, t)/Ot = 0, and Eq. (4) can be written as: 


a Ge 


L Oz Oz 


) ACpp(z) = 0. (10) 


When z —> 0, Cppo is the 7!°Pb,, concentration in the 
interface soils between the plough and the plough pan layer, 
which represents the 7!°Pb,,. concentration at any depth of the 
plough layer or the 7!°Pb,,. concentration in the top soils of the 
plough pan layer. Based on this boundary condition, Eq. (10) 
can be solved as follows: 


Calz) = Cpe OP)? 2, (11) 


Integrating the two sides of Eq. (11) over the depth (z), the 
following can be derived: 


Age = C (D/A). (12) 


By combining Eqs. (5), (6), (11) and (12), the 210Pp con- 
centration in the soils of the plough layer in the steady state 
can be expressed as: 


(13) 


aan ï 
” Ei DEX, 


And the variation of the 7!°Pb,, concentration with depth in 
the plough pan layer in the steady state can be expressed as: 


A i 
C s o TODA 14 
pp (2) H DIA (14) 


By combining Eqs. (9), (12) and (13), the time elapsed 
from initial cultivation to the steady state can be calculated 
as: 


E steady = 1/(An) = 10.37 yr. (15) 


And the diffusion proportion in the steady state can be ex- 
pressed as: 


b=“ x 100% 


(16) 
eN IN, 
H +/D/A 
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where k is the diffusion proportion in the steady state, which 
can be obtained by experiments. 

From Eq. (16), the effective diffusion coefficient of 7!°Pb,, 
is: 


2 
D=\ (=) (17) 


Equations (13) and (14) can be used to simulate the By. k 
depth distribution in the steady state in the soil profile for the 
two layers of the cultivated land. A basic assumption of this 
theoretical derivation is that no erosion or deposition occurs. 
For an eroding soil, the *!°Pb,, inventory would be invariable 
assuming the erosion rate is constant for over 100 years. Then 
the 7!°Pb., depth distribution will also be in steady state, and 
Eqs. (10)-(17) can be applicable as usual. 

True, 7!°Pb,, downward diffusion below the plough depth 
might be partially offset by soil erosion at the surface, caus- 
ing the base of the plough layer to be displaced downwards, 
and some *!°Pb,. that has diffused downwards would return 
to the plough layer. However, the 7!°Pb., distributed below 
the plough depth, more or less, is a matter of fact. In the 
cultivated purple soil of the Sichuan Hilly Basin in China 
(Fig. 1(d)), for example, ca. 37% of the downwards diffused 
210Pb.., according to Eq. (8), returned to the plough layer at 
the first year of diffusion, assuming 7!°Pb.,. is initially quite 
evenly distributed within the plough depth and the soil ero- 
sion rate is 0.5 cm/yr (the real erosion rates are usually less 
than this value [17]). In a word, 210Pp diffusion in cultivated 
soils is a matter of fact in spite of the usual erosion. And the 
proportion would decrease with year and eventually tend to 
constant. Then, the 7!°Pb,, distribution would achieve the 
steady state with a *!°Pb,, tail in the profile below the plough 
depth. 


D. The simulated 7!°Pb,, diffusion process in a soil core from 
the Yimeng mountains area 


To know the 7!°Pb,, depth distribution in cultivated soil 
in reality, a typical 7!°Pb,, profile for a core from the culti- 
vated cinnamon soils of the Yimeng mountains area in China 
was introduced in this study [27], based on which the 7!°Pb,, 
diffusion process could be simulated for the core. The Yi- 
meng Mountains area (ca. 17 180 km?) is located in the cen- 
tral south of Shandong Province, China (34° 22’ to 36°13’N, 
117°24’ to 119°11’E), with elevations within the range of 
150-1165m. The area belongs to the continental monsoon 
climate area of the warm temperate zone, and the cinnamon 
soils developed from limestone weathering crusts are wildly 
distributed. The core was collected from the fields situated 
in flat areas at the top of ploughed fields, where the distur- 
bance except for cultivation over the past 100 years can be 
negligible. The profile (Fig. 2(a)) reflects the steady state of 
210Pp „ depth distribution because the soils have not been dis- 
turbed except for cultivation in nearly ca. 100 years. From 
the profile and the field surveys it can be determined that, the 
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plough depth (H) is ca. 20 cm because the 7!°Pb,,. concentra- 
tions in soils of the upper ~ 20 cm are relatively uniform, sug- 
gesting the *!Pb,, is mixed by tillage within the uppermost 
depth of ~ 20cm, which is also in agreement with our field 
surveys; the 7!°Pb,, diffusion proportion (k) is 32%, which 
is the ratio of the 7!°Pb,, inventory below the plough depth 
to the total inventory in the profile; the effective diffusion 
coefficient (D) is 2.72 cm*/yr, which can be calculated using 
Eq. (17). 

Equations (7) and (8) can simulate the 210Pp redistribu- 
tion process in the plough and plough pan layer of cinnamon 
soils after the initial cultivation before the steady state (e.g. 
t = 2, 4, 6, and 8 yr); Eqs. (13) and (14) can simulate the 
stable *!°Pb.,, depth distribution in the two layers of the culti- 
vated cinnamon soils (t > 10.37 yr) (Fig. 2(b)). 


HI. MASS BALANCE MODEL IN CONSIDERATION OF 
THE DIFFUSION 


A. Model derivation 


Considering *!°Pb,, depth distribution in the cultivated 
land in view of the radionuclide diffusion from the plough 
to the plough pan layer, Eq. (1) can be rewritten as: 


I—\A—AHA,/H =0. (18) 


By combining Eqs. (2), (6), (13), (17) and (18), the erosion 
rate of the cultivated soil can be expressed as Eq. (19) based 
on the relationship between the *!°Pb,, deposition flux and 
the reference inventory (J = \Aver). 


AHI 


Therefore, considering the *!°Pb,. diffusion from the 
plough layer to the plough pan layer, the erosion rates can 
be estimated from Eq. (19). By ignoring the influence 7!°Pb,, 
diffusion (k = 0), Eq. (19) can be reduced to Eq. (3). 


B. Comparison of the revision and existing models 


As above, the difference between the revision and exist- 
ing model lies in whether to take into account the amount of 
210Ph , distributed below the plough depth. To compare the 
two models, let AH; and AH; be the erosion rates estimated 
using Eqs. (3) and (19), respectively, and we have: 


AH, = \HI/(1— J), (20) 


AH» = \H1/[(1—k)(1 — D). (21) 


It is clear that AH, < AA, and the difference ratio (d) of 
AH; and AH; is: 


d= (AH — AH,)/AHb, (22) 
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Fig. 2. (a) The typical 210Db.. profile for a core from the cultivated cinnamon soils of the Yimeng mountains area [27]; (b) The simulated 


210Db, y redistribution process in the core after the initial cultivation. 


By combining Eqs. (20)-(22), one has: 


d=k. (23) 
So, the difference ratio can be represented by the 7!°Pb,, 
diffusion proportion in the steady state of the 7!°Pb., depth 
distribution, and the difference increases with the proportion 
of diffusion. 

If the *!Pb,, diffusion proportion (k) and loss proportion 
(1) are known, Eq. (19) can be used to estimate the erosion 
rates. These parameters can be derived from the experiments. 
Because the k is determined by the 7!°Pb,, profile measure- 
ment, the accuracy of k is related to the accuracy of experi- 
ment, e.g. the study site selection, sampling method, the sam- 
pling section increments, laboratory measurement and etc.. 
Moreover, the recognition of the plough depth, according to 
the shape of the profile in conjunction with field survey, is 
also important to the accuracy of k. 


IV. DISCUSSION AND CONCLUSION 


As above, the erosion rates estimated using the revision 
model are higher than those estimated using the existing 
model by comparison. In other words, the existing model 
underestimates the erosion rate of the cultivated soils to some 
extent because of ignoring the amount of 7!°Pb,, distributed 
below the plough depth. Nevertheless, the existing model in- 
volves the superiorities including compact form, fewer pa- 
rameters, and convenience for application. In contrast, the re- 
vision model involves the limitations including complicated 


form and more parameters for application. Eq. (23) sug- 
gests that the underestimation of soil erosion from the exist- 
ing model is as well less than 10% when the *!°Pb,, diffusion 
proportion (k) is lower than 10%, which could be negligi- 
ble. In other words, the existing model is as usual applica- 
ble to soil erosion estimation when the diffusion proportion 
is lower than 10%. When the diffusion proportion ranges be- 
tween 10%-15%, the existing model is also suitable if the 
underestimation is acceptable in some practical applications. 
However, the underestimation could not be neglected when 
the diffusion proportion is greater than 15%, and the revision 
model would be more applicable to erosion estimation than 
the existing model. No matter what the diffusion proportion 
is, the revision model is universal because the existing model 
is merely a case of the revision model (when k = 0). Ina 
word, the application scope of revision model is universal, 
especially for the locations where the diffusion proportions in 
the cores are greater than 15%. 


Equation (17) indicates that the diffusion proportion is re- 
lated to the diffusion coefficient (D) and the plough depth 
(H), viz. the diffusion proportion increases with the diffusion 
coefficient while decreases with the plough depth. The 
diffusion coefficients and the plough depths vary in different 
areas with different soil types and tillage systems. In the re- 
ported works, the diffusion proportions are usually less than 
10%, e.g. a typical range of 5%-7% at the study sites in 
UK [16, 19] and southern Italy [28, 29]. In these cases, the 
diffusion proportion can be negligible, and the existing model 
is applicable as before. For cultivated brown and cinnamon 
soils of the Yimeng mountains area in China, however, the 
diffusion proportions are greater than 15% [27], so the revi- 
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sion model is more applicable to erosion estimation by com- 
parison. 

Neglecting the 7!°Pb,, below the plough depth, as above, 
will always underestimate soil erosion, nevertheless, it can 
also lead to an overestimation in some cases, e.g., erosion 
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